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Abstract 

The phase diagram of the system OF-LUF3 has been revised using thermal analysis. Specific heat capacity and en- 
thalpy of phase transition and fusion were measured by differential scanning calorimetry for all compounds belonging 
to the system. A thermodynamic optimization of the UF-LUF3 phase diagram was performed by fitting the Gibbs en- 
ergy functions to experimental data that were taken from the literature or measured in this work. Excess energy terms, 
which describe the effect of interaction between the two fluoride compounds in the liquid solution, were expressed by 
the Redlich-Kister polynomial function. The theoretical phase diagram assessed was in suitable agreement with the 
re-evaluated experimental data. 
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1. Introduction 

LiLuF4 (LLF) crystals have been largely investigated 
for doping with rare earth ions, mainly focusing on 
the development of laser media [Ql nana. Taking 
into account that LLF has the most compact crystalline 
structure among the LiREF4 (RE= rare earth element or 
Y, respectively) crystals family, some studies have also 
reported the Li(Lu,RE)F4 (RE = Gd or Y) mixed crys- 
tals as an alternative to achieve better optical properties, 
either wider emission bandwidthsjaa or improvement 
of their photochemical stability |8fl for rare earth doped 
crystals. 

Conversely, there is a lack of information concern- 
ing the thermodynamic properties of LLF and the cor- 
responding LiF-LuF^ phase diagram. Heat capacity 
and enthalpy data for LUF3 were initially investigated 
by Spedding |9] and, more recently, Lyapunov et al. re- 
vised the enthalpy data andproposed a new Cp(T) func- 
tion for this compound 1 1 Oil - No thermodynamic data 
are available for LLF whilst LiF data are given by the 
Barin compilation 111 311 . 

The phase diagram of the binary system LiF-LuF3 
was initially studied in the 1960s by Thoma et al. Ill ill . 



In that work, the LiF-LuF^ phase diagram was de- 
scribed with an intermediary compound (LLF) which 
melts congruently at 1098 K. Two eutectic points were 
settled, one at 22 mol% LUF3 and 968 K and the other 
at 54mol% LUF3 and 1083 K. The polymorphic transi- 
tion from orthorhombic to hexagonal structure in L11F3 
was reported at 1218 K. Harris et al. Il2ll revised this 
phase diagram in the 1980s, determining the melting 
point of LLF at 1123K, the first eutectic was reported 
at 20mol% LUF3 and 977 K and the second one at 
58 mol% LuF 3 and 1 105 K. Both of them used thermal 
analysis to determine the invariant reactions and respec- 
tive temperatures. In this work, the LiF-LuF3 binary 
system was experimentally revised through DSC tech- 
nique. The enthalpy of phase transition and fusion for 
all compounds which belong to this system were cal- 
culated. Using the obtained data and those collected 
from literature IllOL Il3ll - LiF-LuF3 theoretical phase di- 
agram has been optimized. Gibbs excess energy terms 
for the liquid solution, which describe the effects of in- 
teraction between the two fluorides, were expressed by 
the Redlich-Kister polynomial function iflUl . Enthalpy 
of formation and entropy at 298. 15 K were assessed for 
LUF3 and LLF. 
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2. Experimental 

Samples with the desired proportions were mixed us- 
ing commercial LiF (AC Materials, 99.999%) and LuF 3 
(AC Materials, 99.99%), and afterwards were pulver- 
ized in a mortar. DSC curves were obtained using a 
Netzsch STA 409 PC Luxx heat-flux differential scan- 
ning calorimeter. The sample carrier was calibrated for 
T and sensitivity at the phase transformation points of 
BaCC<3 and at the melting temperatures of In, Zn and 
Au. The experiments were carried out under Ar flow of 
30cm 3 /min, using Pt/Au crucibles with lid. Two heat- 
ing/cooling cycles with heating rate of lOK/min and 
samples masses around 30 mg were adopted. The melt- 
ing point of compounds, the temperature of the solid 
phase transition and of the invariant reactions were cal- 
culated considering the extrapolated onset of the ther- 
mal event. The liquidus temperatures of the intermedi- 
ary compositions were evaluated from the extrapolated 
offset temperatures. 

The heat capacity measurements were performed 
using the same equipment described above, where a 
proper sample carrier for Cp analysis was installed. The 
experiments consisted of three steps, an isothermal seg- 
ment at 40 °C for 20min, a dynamic heating segment 
with heating rate of 10 K/min and a final isothermal seg- 
ment at maximum temperature for 5 min. Other exper- 
imental conditions were kept as mentioned previously 
for simple DSC experiment. LuF 3 and LiLuF4 heat ca- 
pacity were determined according to ratio method. In 
this method, the sample DSC heat flow signal is com- 
pared to the DSC signal of a calibration standard of 
known specific heat (sapphire in this case). Both curves 
are corrected by a baseline correction experiment where 
empty reference and sample crucibles are placed in the 
DSC furnace and the system signal drift is measured un- 
der identical experimental conditions. 

3. Thermodynamic method 

A thermodynamic simulation of a T - X binary phase 
diagram requires the description of Gibbs energy func- 
tions for all compounds in the system and the Gibbs 
functions of mixing, if solution phases are present. In 
most cases the excess Gibbs functions are unknown for 
the solutions, therefore a thermodynamic assessment is 
necessary in order to determine thermal effects of mix- 
ing. Gibbs energy for a solid compound is defined as a 
function of enthalpy and entropy at the reference tem- 
perature state (298.15 K) and can be obtained from the 
heat capacity function (Cp(T)) as follows: 
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Combining the equations (f2]i and ©, the G(T) equa- 
tion is given as follows: 
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Usually the heat capacity is the most accessible phys- 
ical property to be measured. Cp(T) functions can be 
determined by fitting a set of experimental data at a 
proper polynomial function (e.g. equationQ]). Depend- 
ing on the compound, more terms are added or disre- 
garded in the polynomial function to obtain the best 
fitting. According to equation (Q} and considering the 
equations © and (O, the enthalpy of formation, the ab- 
solute entropy at the reference temperature and the heat 
capacity are required to determine the minimum of G 
and thus the thermodynamic equilibrium. These data 
are not available for the LLF, therefore the experimental 
data calculated by DSC was assumed to set the Cp equa- 
tion, as will be discussed later. Cp equation for LuF 3 
was also obtained using this technique and was com- 
pared with the literature values. Enthalpy and entropy at 
298. 15 K for LLF and LuF 3 were properly assessed by 
optimization. Table [T] summarizes all calorimetric data 
used on the phase diagram assessment. The thermody- 
namic data for the LiF end member is well established 
and was taken from Barin compilation 111 311 . 

For the liquid solution phase, function G(T) is ex- 
pressed as the sum of the Gibbs energy weighed con- 
tribution of the pure compounds (Go), the contribution 
of an ideal mixture (Gjd) and finally a term related to 
the non-ideal interaction, defined as the excess energy 
(G ex ). The sub-regular solution model of Redlich-Kister 
was adopted to describe the excess energy of the liquid 
phase in this system fl4Tl . and is given by: 

N 

G ex = x A x B 2^ Lj (x A - x B ) J (5) 
J=o 

where x A and x B are the molar fractions of components 
A and B, respectively. Lj terms represent the interac- 
tion coefficients between the basis compounds and they 
are given as a linear function of temperature. The opti- 
mization was performed using the OptiSage module in 
the FactSage 6.2 software [15], which uses the Bayesian 
Algorithm lfl6ll . This algorithm is based on a probabil- 
ity model to obtain the fit between the theoretical Gibbs 
energy functions and the experimental data. 
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Table 1 : A#(298. 15 K) (kJ mor 1 ), S (298. 15 K) (J K _1 mor 1 ), AH PT (Phase transition: heat of fusion or polymorphic transition) (kj mor 1 ) and 
Cp data a, b, c entering eq. {T] for LiF, L11F3 and the intermediate compound LiLuF4. 



Compound 


A//(298.15/Q 


5(298.15/:) 


AHp T (Lit) 


AHpj a 


a 




b 


C 


LiF(S)^ 


-616.931 


35.660 


27.09 


27.68 


42.689 


1.742 x 10 


-2 


-5.301 x 10 5 


LiF(l) b 


-594.581 


42.997 






64.183 








LuF 3 (Sl) fl 


-1697.662 


118.420 


25.40^ 


18.28 


97.496 


9.460 x 10 


-3 


-9.998 x 10 5 


LuF 3 (S2) fl 


-1693.056 


111.358 


29.90 c 


23.34 


114.500 c 








LuF 3 (/) fl 


-594.581 


42.997 






131.800 e 








LiLuF^j 


-2323.639 


153.218 




63.538 


139.430 


3.200 x 10 


-2 


-1.881 x 10 6 



"AHpT and Cp function were measured by the DSC technique; A//(298.15 K) and 5(298.15 K) were assessed in this work. 

*Data taken from Barin H3l . 

r Data taken from Lyapunov et al. fl9l . 



4. Results and discussion 

4.1. Experimental results 

DSC curves for pure LuF 3 and six different composi- 
tions are shown in Figure Q] The curves for 17, 20 and 
27 mol% LuF 3 exhibit one endothermic peak with onset 
near to 968 K. This peak is due to the L ^ LiF+LiLuF4 
eutectic, and the 20mol% LuF 3 DSC curve shows a sin- 
gle peak at this temperature, since the eutectic point is 
very close to this composition. The second broader en- 
dothermic peaks at 17 and 27mol%LuF 3 curves mark 
the end of melting of the primary phase (LiF and LiLuF 4 
respectively) at the liquidus. 

Curves with compositions of 50 and 58 mol% il- 
lustrate the LLF congruent melting at 1126K, and the 
second eutectic reaction at 1115K, respectively. DSC 
curve for 80mol%LuF 3 shows, beyond the small peak 
representing the liquidus line, also two intermediate 
peaks with onset at 1115 K and 1216 K corresponding 
to the second eutectic reaction and the LuF 3 polymor- 
phic transition from orthorrombic to hexagonal, respec- 
tively. 

Finally, in the LuF 3 curve the polymorphic transition 
(1216 K) and fusion (1403 K) can be recognized very 
clearly. The main features of the LiF-LuF 3 experimen- 
tal phase diagram are summarized in Table |2] together 
with some previously reported data. In general the data 
obtained in this work are in agreement with those from 
Harris et al. phase diagram 11211 . 

The heat of fusion and polymorphic transition 
(AHpj) for LiF, LuF 3 and LiLuF 4 were calculated from 
the DSC peak area of these thermal effects for each 
compound. In a previous work, the experimental AHpr 
for LiF was reported and reasonable agreement with the 
literature value was found 11131 Il7ll . Compared to the lit- 
erature data, the main difference in the results were no- 
ticed for LuF 3 thermal events. The evaluated AHpr for 
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Figure 1: Second heating DSC curves for some compositions in the 
LiF-LuF} system. 



LUF3 in this work are approximately 20 % smaller than 
those reported by Lyapunov et al. ft 1 Ofl . Since our DSC 
experiments were repeatedly confirmed and the raw ma- 
terials used were very pure, these new results have been 
taken into account to the final calculation of the LiF- 
LUF3 phase diagram. The heat of fusion and heat ca- 
pacity (Cp function) for the LiLuF4 compound have not 
yet been reported in the literature. Therefore, for the 
first time these calorimetric properties were evaluated. 
Figure[2]compares the Cp data calculated through DSC 
technique and those estimated using the corresponding 
Cp(T) functions of end members compounds and con- 
sidering the Neumann-Kopp rule. One can see, for the 
considered temperature range, there are no signicant de- 
viations (not bigger than 3%) between the experimental 
and estimated Cp data. Furthermore, taking into ac- 
count the polynomial fitting of the DSC experimental 
data using the function (Q]i, the a, b and c parameters 
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Table 2: Main features of LiF- 


-L11F3 experimental phase 


diagram 


and the previously reported data. 






Composition (Mol 9 


LuF 3 ) 




Temperature (K) 




Equilibrium Reaction 


This work 


Harris 


12} Thoma [11J 


This work 


Harris [JJJ 


Thoma [jTJ 


L ^ LiF+LLF 


20 


20 




22 


968 


977 


968 


(Eutectic) 
L^LLF 


50 


50 




50 


1126 


1123 


1098 


(Congruent melting) 
L ^ LLF+LuF 3 (5 1) 


58 


58 




54 


1115 


1105 


1083 


(Eutectic) 

L+LuF 3 (5 1) f± L+LuF 3 (52) 










1216 




1218 


(Polymorphic transition) 

















obtained define the Cp(T) function for the OLUF4 com- 
pound. 

Heat capacity re-evaluation for the L.UF3 orthorhom- 
bic phase was also carried out by DSC (Figure [3}- The 
results are in reasonable agreement with those reported 
previously by Spedding et al. Some small de- 

viations from the general behavior of curve were ob- 
served for the measured Cp data (e.g. 930 and 1065 K). 
Those undesirable effects are due small fluctuations on 
the DSC base line and they can be disregarded. Taking 
that aside, the difference between the evaluated and re- 
ported data is not larger than 4%, even for the higher 
temperatures considered. Therefore, the Cp(T) param- 
eters set by fitting a polynomial model ([TJ on DSC ex- 
perimental data were assumed to perform the numerical 
optimization in the LiF-LuF3 system. 
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C p model: C p (T) = a + b7 + c/T 
Chi A 2/DoF = 0.24377 
R A 2= 0.99744 
a = 139.4 ±0.4 
b = (32 ± l)E-3 
c=(-188±3)E+4 
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Figure 2: LiLuF4 Cp data measured experimentally through DSC 
technique and Cp function estimated by Neumann-Kopp rale. 



4.2. Thermodynamic assessment 

Once the T —X experimental points are available and 
the calorimetric data for each compound belonging to 
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Figure 3: L11F3 Cp data obtained experimentally through DSC tech- 
nique and Cp data reported previously by Sppeding 3. 



the system were evaluated or collected from the liter- 
ature, a proper optimization of the binary LiF-LuF3 
phase diagram became possible. The excess Gibbs 
energy for liquid phase has been optimized accord- 
ing to the Redlich-Kister polynomial model using the 
Bayesian Optimization Algorithm of FactSage 111511 . 

Figure |4] presents the LiF-LuF3 optimized phase di- 
agram together with the T - X experimental points for 
comparison. It may be highlighted an excellent agree- 
ment of the assessed temperature values for the eutec- 
tic reactions, polymorphic transition and even for the 
liquidus line. L1LUF4 melting point was confirmed at 
~1 123 K. The value calculated for LiF/LiLuF4 eutectic 
composition was about 3 mol% larger than the experi- 
mental value. However, it should be noticed that, exper- 
imentally is pretty hard to separate close thermal events, 
as it happens with liquidus and solidus line nearby to 
the eutectic composition. In general, the experimen- 
tal determination of the liquidus from offsets is much 
less accurate than the determination of eutectic or phase 
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transition temperatures. The assessed enthalpy and en- 
tropy at 298. 15 K temperature for the LuF^ (orthorhom- 
bic and hexagonal phases) and for LiLuF4, are listed in 
table Q] The excess parameters calculated, given by Lq 
and L\ in the Redlich-Kister model (equation [5]) were 
L = -69089.24 + 49.37F and U = 8903.63 - 1.587. 
Even though those numbers can not be taken as unique, 
they represent the best assessment achieved. Together 
with the other data presented in this work, that numer- 
ical simulation gives a complete thermodynamic de- 
scription for the LiF-LuF3 binary phase diagram. 
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Figure 4: Assessed phase diagram of the system L1F-L11F3 (lines) 
together with the experimental data from DSC (squares). 



5. Conclusions 



Thermodynamic assessment has been performed on 
the LiF-LuF3 binary system and the excess Gibbs en- 
ergy terms for the liquid phase could be properly op- 
timized using the Redlich-Kister polynomial model. 
Based on our experimental results and the thermody- 
namic assessment it was possible to confirm the values 
for the two eutectic temperatures within this system. A 
Cp(T) function of L1L.UF4 compound has been set for 
the first time and it was in accord with Neumann-Kopp 
rule. It has been successfully confirmed the previously 
Cp data found in the literature for LUF3 (orthorhom- 
bic phase). A//(298.15K) and S (298. 15 K) have been 
assessed for L.UF3 and LiLuF4, and a re-evaluation for 
AHpp data in the LUF3 compound was performed. This 
work offers a more complete thermodynamic descrip- 
tion of L1F-LUF3 system and makes it possible further 
thermodynamic assessment for ternary and multi com- 
ponent phase diagram based on this binary one. 
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